MODULE CBALANCE_CONST
INCLUDE 'soil.inc'
     real(8), PARAMETER :: a15(nv2) = (/4.0, 2.53, 2.77, 48.0, 48.0, 9.0, 100.0 /)
      real(8), parameter :: ab_litter(nv2) = (/0.59, 0.58, 0.6, 0.64, 0.295, 0.08, 0.71 /)
      real(8), parameter :: n_littera(nv2) = (/0.91, 0.76, 1.4, 1.5, 1.25, 2.5, 0.82/)
      real(8), parameter ::  ash_littera(nv2) = (/2.13, 1.5, 4.6, 3.3, 4.3, 8.2, 1.4 /)
      real(8), parameter ::  n_litterb(nv2) = (/0.71, 0.56, 0.6, 0.64, 0.63, 1.3, 0.76/)
      real(8), parameter :: ash_litterb(nv2) = (/1.0, 0.75, 2.4, 1.6, 0.64, 3.0, 0.89 /)
! monthly leaf area index + stem area index, one-sided
! monthly leaf area index, one-sided
      real(8) tai(nv2,12), gai(nv2,12)

 

!top of canopy (m)
      real(8), parameter :: htop(nv2) = (/35.0,20.0,17.0, 0.5,&  
     &             0.5, 0.5, 0.5/) !, 0.0/

!bottom of canopy (m)
      real(8), parameter :: hbot(nv2) = (/ 1.00,11.50, 8.50, 0.01,&  
     &             0.10, 0.01, 0.01/) !, 0.00/
! minimum temperature for photosynthesis (kelvin)
      real(8), parameter :: tmin(nv2) = (/278.16,273.16,268.16,273.16,&
     &           273.16,273.16,273.16/) !,  0.00/
 
! saturated soil albedos: 1=vis, 2=nir
!      real(8), parameter ::(albsat(i,1),i=1,9)/0.12,0.11,0.10,0.09,0.08,0.07,0.06,0.05,0.22/
!      real(8), parameter ::(albsat(i,2),i=1,9)/0.24,0.22,0.20,0.18,0.16,0.14,0.12,0.10,0.34/

! dry soil albedos: 1=vis, 2=nir 
!      real(8), parameter ::(albdry(i,1),i=1,9)/0.24,0.22,0.20,0.18,0.16,0.14,0.12,0.10, 0.34/
!      real(8), parameter ::(albdry(i,2),i=1,9)/0.48,0.44,0.40,0.36,0.32,0.28,0.24,0.20, 0.58/

! leaf reflectance: 1=vis, 2=nir
      real(8)  rhol(2,nv2), rhos(2,nv2),taul(2,nv2),taus(2,nv2)
     real(8), parameter :: betads = 0.5, betais = 0.5

! leaf/stem orientation index: valid range = -0.4 to 0.6
      real(8), parameter :: xl(nv2) = (/ 0.10, 0.25, 0.01,-0.30, & 
     &          0.25, -0.30,-0.30/) !, 0.00/
     ! minimum temperature for photosynthesis (kelvin)

!photosynthetic pathway: c3 = 1, c4 = 0
      real(8), parameter :: c3psn(nv2) = (/1., 1., 1., 1., 1., 1., 1./) !, 1./

! co2 michaelis-menten constant at 25c (pa)
      real(8), parameter :: kc25(nv2) = (/30., 30., 30., 30., 30., 30., 30./)

!q10 for kc25
       real(8) akc(nv2)
! o2 michaelis-menten constant at 25c (pa)
      real(8), parameter :: ko25(nv2) = (/30000.,30000.,30000.,30000.,30000.,30000.,30000./)

! q10 for ko25
      real(8) ako(nv2) 

! maximum rate of carboxylation at 25c (umol co2/m**2/s)
      real(8), parameter :: vcmx25(nv2) = (/50.,33.,33.,33.,17.,&
     &             33.,50./) !, 0./
!       real(8), parameter :: vcmx25 /7*250./
! q10 for vcmx25
      real(8) avcmx(nv2)  
! minimum leaf conductance (umol/m**2/s)
      real(8), parameter :: bp(nv2) = (/8210.,8210.,8210.,8210.,8210.,8210.,8210./) ! 7*2000./

! slope for conductance-to-photosynthesis relationship
       real(8), parameter :: mp(nv2) = (/9.,9.,6.,9.,9.,9.,9./) !,9./
!      real(8), parameter :: mp /7*9.0/ 
! quantum efficiency at 25c (umol co2 / umol photon)
      real(8), parameter :: qe25(nv2) = (/0.06,0.06,0.06,0.06,0.06,0.06,0.06/) !, 0.00/

! q10 for qe25
      real(8), parameter :: aqe(nv2) = (/1.0,1.,1.,1.,1.,1.,1./)
!foliage nitrogen concentration when f(n)=1 (%)
      real(8), parameter :: folnmx(nv2) = (/1.5,1.5,1.5,1.5,1.5,1.5,1.5/) !, 0.0/
! empirical parameter that defines root fraction decrease with depth
      real(8), parameter :: rdp(nv2) = (/0.94,0.94,0.94,0.97,&
     &          0.97,0.94,0.94/) !,1.00/
!top of canopy (m)

      real(8), parameter :: dleaf(nv2) = (/0.04,0.04,0.04,0.04,0.04,0.04,0.04/)
      real(8), parameter :: z0m(nv2) = (/2.62, 1.1, 0.94,0.06, &
          & 0.06,0.06,0.06/)
          real(8), parameter :: dd(nv2) = (/23.45, 13.4, 11.39, 0.34, &
          & 0.34, 0.34, 0.34/)
          real(8), parameter :: r_a(nv2) = (/7.0, 6.0, 7.0, 11.0, 7.0, 7.0, 6.0 /)
          real(8), parameter :: r_b(nv2) = (/1.0, 2.0, 2.0, 2.0, 1.5, 1.5, 3.0 /)

!--------------------------------------------------------------------

END MODULE CBALANCE_CONST
